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The complex iron line in NGC 7469 observed by BeppoSAX 
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ABSTRACT 

In this letter we present analysis of BeppoSAX data from a long look at the Seyfert 1 galaxy 
NGC 7469. The presence of a soft excess below 0.8 keV is confirmed by our analysis and no 
warm absorber component is required to fit the spectrum. A complex iron emission line and 
the Compton reflection hump are clearly detected. The profile of the line is too broad to asso- 
ciate this feature with distant matter. In addition, the observed soft excess and the energy of the 
iron line Ei? e =6.8 keV strongly support a scenario in which the hard X-rays are reprocessed 
by a photoionized accretion disc. This hypothesis was tested fitting the BeppoSAX spectrum 
with the ionized disc reflection model of Ross & Fabian. A second narrow line component, 
in addition to that produced in the disc, is also required to fit the observed iron line profile. A 
high energy cut-off around 150 keV is clearly detected in the spectrum. 
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1 INTRODUCTION 



A narrow iron line component is clearly detected in recent Chan- 
dra (Kaspi et al. 2001, Yaqoob et al. 2001) and XMM observations 
(Reeves et al. 2001, Pounds et al. 2001) of some Seyfert 1 galax- 
ies. These observations suggest that part of the Fe line is produced 
very far from the accretion disc. A strong broad iron line compo- 
nent is found in the X-ray spectra of many Seyfert galaxies (Fabian 
et al. 2000) by ASCA (Nandra et al. 1997, Tanaka et al. 1995, 
Yaqoob et al. 2002), XMM (Wilms et al. 2001) and BeppoSAX 
(Guainazzi et al. 1999, De Rosa et al. 2002). 

The near (z~0.017) Seyfert 1 NGC 7469 shows a very com- 
plex X-ray spectrum. EXOSAT observed a soft excess in this source 
(Barr 1986) which was confirmed by Einstein (Turner et al. 1991). 
Two ROSAT observations in 1991 (Turner et al. 1993) and 1992 
(Brandt et al. 1993) did not resolve the origin of the soft spec- 
tral component (e.g. a soft excess, a warm absorber or a combi- 
nation of the two). No evidence of a warm absorber was found in 
the 1993 ASCA spectrum (Guainazzi et al. 1994). An iron emis- 
sion line and a flattening of the spectrum above 10 keV were de- 
tected by Ginga (Piro et al. 1990). When the spectrum above 10 
keV was reproduced with a cold Compton reflection component 
(George & Fabian 1991; Matt, Perola & Piro 1991), the intrinsic 
photon index was Y= 1.99 zztg; 8s (Nandra 1991), while a partial 
covering model yielded Y= 1.92 ± 0.03 (Piro et al. 1990). The 
line profile was found to be narrow (a < 150 eV) in the ASCA 
spectrum (Guainazzi et al. 1994), suggesting that the reprocess- 
ing material has to be far from the central source (several hun- 
dreds of Schwarzschild radii). Nevertheless Nandra et al. (1997) 
found marginal evidence of a broad component of the iron line in 
analysing the ASCA spectrum. 

The broad band of BeppoSAX (Boella et al. 1997) is partic- 
ularly suited to the deconvolution of such complex spectra. In this 



letter we present the analysis of a BeppoSAX week-long observa- 
tion of NGC 7469 in 1999. This spectrum was already analysed by 
Perola et al. (2002) within a sample of nine sources observed with 
BeppoSAX. Their analysis showed evidence of soft excess com- 
ponent and a not resolved iron line was also detected, indicating a 
component from the centre of the accretion flow. 



2 OBSERVATION 

BeppoSAX observed NGC 7469 from 1999 November 23 to 1999 
November 29. The net exposure time in the MECS was 249180 s. 
The observed mean flux was F2_iofc e v = 3.7 x 10 _11 erg cm~ 2 
s _1 for our best fit model (see Table [l]). Spectra were extracted 
within a circular region centered on the source with radii of 4' 
and 6' for MECS and LECS respectively. The background was ex- 
tracted from event files of source-free regions ("blank fields"). The 
PDS spectrum was filtered with fixed rise time. The BeppoSAX 
data were then fitted using the XSPEC 11.1 package. All quoted 
uncertainties correspond to 90 per cent confidence intervals for 
one interesting parameter (Ax 2 =2.71). Each model we tested was 
multiplied by a normalization constant in order to take into ac- 
count possible miscalibrations between the different instruments. 
We allowed the PDS/MECS normalization to vary between 0.77 
and 0.95, while the LECS/MECS normalization ratio was running 
between 0.7 and 1 (Fiore, Guainazzi & Grandi 1999). 

In Figure [l] we show the LECS (0.7-2.5 keV), MECS (5-10 
keV) and PDS (13-200 keV) lightcurves, the hardness ratio MECS 
(5-10 keV) / MECS (3-10 keV) is also plotted in the last panel. In 
each plot we reported the x 2 / dof for a constant hypothesis. Even if 
the source showed flux variability in the soft, medium and hard en- 
ergy range (P x 2 < 10" 3 for LECS, MECS and PDS lightcurves), 
the MECS hardness ratio (which is in the energy range of interest 
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BeppoSAX. NGC 7469 (1999) 



NGC 7469. Iron emission line 
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Figure 1. LECS. MECS, PDS lightcurves and MECS hardness ratio of Bep- 
poSAX observation of NGC 7469. In each panel is indicated the x 2 /dof for 
a constant hypothesis. 

NGC 7469. Absorbed power law mode 
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Figure 2. LECS, MECS and PDS data (upper panel), and data/model ratio 
(lower panel) in the case of an emission continuum fitted with a power law 
absorbed by a cold galactic gas. 



to investigate the Fe line behaviour) is in good agreement with a 
constant value (P x 2= 0.33). A detailed discussion about the com- 
plex spectral variability in NGC 7469 (Nandra & Papadakis 2001 
and references therein), is beyond the focus of this letter. 



3 RESULTS 

We fitted the LECS (0.2-3 keV), MECS (1.5-10 keV) and PDS (13- 
200 keV) data simultaneously. In Figure |^ we show the data and 
data/model ratio when the whole BeppoSAX spectrum is fitted by 
an absorbed (Nh = 4.8 x 10 20 crn~ 2 the Galactic value from Elvis 
et al. 1989) power law. This model is clearly inadequate, with sev- 
eral features emerging from the spectra: an excess at energy less 
than 1 keV, an iron line component with the Compton reflection 
hump, and a deficit of counts above 70 keV. 

The baseline model we employed to fit the data consisted of an 
absorbed cut-off power law, a cold reflection component (PEXRAV 
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Figure 3. Ratio data/model of the iron line when the data are fitted with the 
baseline model in the whole BeppoSAX band except the 4-7.5 keV energy 
range. 

NGC 7469. Ionized disc reflection model. Iron line 
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Figure 4. The unfolded BeppoSAX/MECS spectra of NGC 7469 is plotted 
with the ionized disc reflection model with solar Fe abundance (see Table 
bJ). Clear residuals around the Fe line can be observed. 



NGC 7469. Ionized disc reflection model + Gaussian line 
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Figure 5. A similar plot to Figure ^ but with the spectrum fitted with an 
ionized reflection model plus a gaussian component. 
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Table 1. Cold disc reflection model. M The inclination angle is linked to that of the reflection component and frozen to 30°. (•) Gravitational radius 
Rg=GM/c 2 . W In this model for the narrow component E^=6.3 keV and <r=0 keV. All the line energies are quoted in the source rest frame. 



Model 
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(keV) 


(keV) 
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Single gaussian 


80±f 


j 07+0.02 


1.49±0.40 


245±f?f 


<: 42+0.11 
°-^ z=t 0.12 


>0.39 


175+80 


151.3/136 


Diskline b 
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>52 
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155.5/136 


Diskline 1 ' + Narrow gaussian* 


84± 7 


2.06±0.02 


i.8±8:S 


140±|| 5 


6.82+0.15 


> 15 


121 ± 100 
75 + 40 


142.9/134 



Table 2. Ionized disc reflection model. The exponent for the disc emissivity law r' 3 is frozen to the value f}=-3. W Parameter fixed duringthe fit. W Parameter 
pegged at upper limit. 



Model T log ^ R 6 t rin i out E L X 2 /dof 

(°) (R G ) (Rg) (keV) 

Ionized disc 2.01+0.02 1.08±g;J| 0.95±g;g| 20.6±|;g 6 *10 7 169/139 

Ionized disc + gaussian 2.02+ 0.02 1.21+0.02 0.87±g;?g 37.8+f; 7 , 6 < 51 6.30+ 0.05 151.9/137 

ASCA 



Ionized disc + gaussian 2.02 + 0.03 1.08±g ^| 1 (frozen) <40 6 *10 7 6.36+0.06 399.3/378 



in XSPEC Magdziarz & Zdziarski 1995) and a black body to re- 
produce the excess at low energy. The probability to exceeding F 
with the addition of the high energy cut off in the baseline model 
is Pf ~ 96.5 per cent. In Figure H the MECS data/model ratio is 
shown when the baseline model is applied to the whole BeppoSAX 
band, but excluding 4-7.5 keV energy range where a relativistic iron 
line is expected to be present. The ratio in Figure || shows clearly 
the presence of a complex iron line. We tried several models to fit 
the line profile and all the best fit parameters are shown in Table [l| 
A single gaussian is clearly inadequate (x 2 /dof=151. 3/136). The 
second step was to substitute an unphysical broad gaussian com- 
ponent with a physical model for a relativistic profile of the line 
(DISKLINE in XSPEC, Fabian et al. 1989). The model is still in- 
adequate to reproduce the line residuals (x 2 /dof=155. 5/136). The 
addition of a second narrow gaussian component at the energy of 
the neutral iron yields a significant decrease of \ 2 (Ax 2 =13 with 
the addition of 2 free parameters. Pf=99.70 per cent). The EW was 
121 ± 100 eV and 75 ± 40 eV for the broad and narrow compo- 
nents respectively. The centroid of the relativistic line El=6.82 ± 
0.15 keV, indicates He-like iron. 

Two gaussian narrow lines give a fair x 2 (x 2 /d°f=l42/135). 
The second line (EW=57±36 eV) lies to the energy of highly ion- 
ized iron E=6.94±(];i2 keV, i.e. H-like iron. This component has 
to be produced in a gas characterized by an ionization parame- 
ter £ = 4nL ion /n e R 2 > 2000 erg cm s _1 (Matt et al. 1996). 
Assuming a temperature T=10 6 K, the total recombination rate to 
all states for FeXXVI is a to t = 2.4 x 10~ n cm 3 s" 1 (Seaton 

2 

1959), and the line emissivity is Ju ne = j^AFe yhvottot, with 
Afb and y the iron aboundance and the iron effective fluores- 
cent yield respectively. In the case of a simple spherical geome- 
try if we assume a ionizing luminosity Li 0n ~ 10 43 erg s _1 , and 
with the observed H-like iron line luminosity Lu„ e ~ 10 41 erg 



s 1 , we can estimate a lower limit to the electron number density 

ne > ( , 3 t."" e ) 2 x (%n.)~ 3 = 2.4 x 10 10 cm" 3 (as- 

suming Af b = 3.3 x 10~ 5 of Morrison & McCammon 1983 and 
y ~ 0.63 by Matt et al. 1993b for a temperature of 10 6 K) and then 
an upper limit for the distance of the gas from the ionizing central 
source, R < 7 x 10 14 cm. This means that such gas must lie within 
few hundred Schwarzschild radii (R 3 =2GM/c 2 ) of the central black 
hole (assuming a mass for it of ~ 10 7 Mq, Wandel et al. 1999). 

The fit with two gaussian emission lines (one narrow) is 
still perfectly acceptable (x 2 /dof=133/134), nevertheless the ex- 
tremely large value of the intrinsic width of the broad component 
(<t=1.8±J;°. keV) indicates a poor description of the continuum at 
the energy of the line. Thus we are lead to believe that a diskline 
+ gaussian component is a more suitable and physically consistent 
model to reproduce the Fe emission line observed in NGC 7469. 

The presence of an ionized iron line and of a soft emission 
component, together with an unusual value of the relative reflec- 
tion, significant greater than 1 (R=1.8±o;|, see Table [l]), strongly 
suggests the hypothesis that all these features could be the re- 
sult of hard X-ray reprocessing in a hot ionized accretion disc 
(Matt et al. 1993a, 1996, Ross, Fabian & Young 1999, Nayakshin 
et al. 2000, Ballantyne et al. 2001). We tested an ionized reflection 
model which takes into account the reflected spectrum of the disc 
in the whole BeppoSAX energy range. The model we employed 
was that of Ross & Fabian (1993). The most important quantity in 
determining the shape of the reflected continuum is the ionization 
parameter £ = iirFx /uh, where Vx is the X-ray flux (between 
0.1-100 keV) illuminating a slab of gas with solar abundances and 
constant hydrogen number density n_fr=10 15 cm" . The incident 
flux is assumed be a power law with spectral photon index F and 
a high energy cut-off at 100 keV. The computed reflected spectrum 
is multiplied for a factor R ("reflected fraction") and added to the 
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primary continuum. This model includes the Fe Ka emission line 
and the spectral features at low energy (emission lines and recom- 
bination continuum). Larger values of £ indicates a more ionized 
disc which will affect the strength and width of the iron line (Matt 
et al. 1993a, 1996), and the absorption edges. The soft emission 
from the accretion disc is not take into account to determine the 
reprocessing features. During the fit we applied to the spectrum a 
relativistic blurring appropriate for a Schwarzschild geometry as- 
suming a disc emissivity law (Fabian et al. 1989). 

Fitting the spectrum with an ionized disc model (see Table 
^ with solar Fe abundance does not give a satisfactory result 
(X 2 /do/=169/139). Clear residuals can be observed at the energy 
of the iron line (see Figure [J). To increase the Fe emissivity to re- 
produce the observed strength of the line we tried to fit the data with 
an ionized reflection model with 2 and 5 times solar Fe abundances. 
Nevertheless these models are inadequate to fit the soft part of the 
spectrum, and the result of the fit is poor (x 2 /do/( 2 xFe) = 187/139, 
X 2 /do/( 5 xFe) =205/1 39). A better model to reproduce the iron line 
residuals is to add a narrow line component. The ionized disc model 
with an additional gaussian component finally give us a good fit 
(Ax 2 =17 for the sum of two interesting parameters. Pf=99.97 per 
cent) and no large residuals can be found at the energy of the iron 
line (see Figure ra). The best fit parameters for the several ionized 
reflection models we employed are shown in Table ^. In Figure ^ 
we plot the confidence contours (68, 95 and 99 per cent) Reflection 
fraction and photon index for the ionized disc plus gaussian line 
model. 

No additional (cold or ionized) iron edge is required to fit the 
spectrum. If we add an absorption edge with the energy fixed to the 
value for cold iron E e( j ffe =7.1 keV, the fit worsens and we find a 
very low upper limit for the optical depth r < 0.01. 

The soft excess at low energy can not be accounted for by 
a warm absorber component. If we substitute the black body in 
the last model in Table [l], with two absorption edges at the en- 
ergy of the OVII (0.74 keV) and OVIII (0.87 keV), the fit is worse 
(X 2 /do/=150/134). The same result is obtained using a more de- 
tailed description of the warm absorber with ABSORI model in 
XSPEC. A warm absorber component is not required to fit the low 
energy data of NGC 7469 (Pf=34 per cent). 

Piro et al. (1990) showed that the GINGA spectrum of 
NGC 7469 observed in 1988 could be reproduced by a "partial 
covering" model (Holt et al. 1980). In this scenario an absorbing 
medium is supposed to cover only a fraction f cov of the central 
source with a column density N_ff col) . The direct X-ray would be 
blocked at low energy, but would penetrate through the absorbing 
gas above 5 keV imprinting the curvature in the observed spec- 
trum. Nevertheless Leighly et al. (1996) showed that this model 
was not able to explain the observed flux and spectral variability in 
NGC 7469. 

When we model the BeppoSAX spectrum of NGC 7469 with 
partial covering the fit is fair (x 2 /do f =142. 7/136), with Nh co „ = 
(5 ± 1) x 10 23 cm" 2 and f cov = (0.32 ± 0.04), and the EW of 
the iron line EW=69 ± 20 eV. Nevertheless the expected value of 
the equivalent width for an optically thin medium illuminated by 
an isotropic continuum, with a Fe solar abundance and the Nhcov 
and f cov obtained by BeppoSAX observation, is in the range 100- 
200 eV, i.e. larger than that observed in the spectrum. This evidence 
together with the complex profile of the Fe line and the soft X-ray 
excess, make the reflection from a (possibly ionized) accretion disc 
the most appealing scenario. 



4 DISCUSSIONS AND CONCLUSIONS 

The BeppoSAX long look at NGC 7469 shows a very complex 
X-ray spectrum. Both a soft excess component and reprocessing 
features are detected. The continuum emission can be reproduced 
either with a cold reflection model and a black body component 
(with temperature kT ~ 84 eV fully consistent with that measured 
by ASCA, Guainazzi et al. 1994) or with an ionized disc reflection 
model that takes into account the disc emissivity below 0.8 keV. 
Whichever model we employ to fit the spectrum, the observed iron 
line requires a complex profile which is hard to associate with a 
distant reflector. In addition the energy of the relativistic line E=6.8 
keV corresponds to high ionized iron that can not be reconciled 
with the cold disc scenario. 

An ionized iron line component was already observed by 
XMM in Mkn 509 (Pounds et al. 2001), Mkn 205 (Reeves 
et al. 2001) and NGC 5506 (Matt et al. 2001). In this last source a 
solution in term of a blend of He- and H- like narrow iron line is 
preferable to that of a relativistic ionized disc, but in Mkn 509 and 
Mkn 205 this feature was explained in terms of an ionized accre- 
tion disc. In the case of NGC 7469 the predicted iron line compo- 
nent from the ionized disc model does not fit the whole observed 
profile. A narrow line component has to be included in the model. 
The origin of this second component is still a matter of debate. It 
has been detected in many Seyfert Is observed so far by Chan- 
dra (Kaspi et al. 2001, Yaqoob et al. 2001) and XMM (Pounds 
et al. 2001, Reeves et al. 2001). It is associated with a region dis- 
tant from the central source and probably coincident either with a 
Compton thick gas, e.g. the "torus" obscuring the Seyfert 2s (An- 
tonucci 1993, Ghisellini et al. 1994) or with the Compton thin 
BLR. At least in the cases of NGC 4051 (Guainazzi et al. 1998) 
and NGC 5506 (Matt et al. 2001) the narrow core of the iron line 
was associated to a Compton thick gas. Nevertheless the bulk of 
the Compton reflection in NGC 7469 can be completely associated 
with the ionized disc (see Figure []) iron component. In fact with 
the combination of a cold and an ionized Compton reflection we 
obtained a very poor fit (x 2 /do/=156.4/135). 

Even if a partial covering model can reproduce the observed 
BeppoSAX spectrum, the value of the equivalent width of the iron 
line predicted in this model is significantly larger than that mea- 
sured. This could require that the covering medium has a small 
solid angle at the source. 

A warm absorber component is not required to fit the Bep- 
poSAX low energy spectrum, this allowed us to detect the redward 
extent of the iron line, and thus constrain the inner radius of the 
disc. 

A double iron line component also can reproduce the observed 
profile. We cannot rule out the possibility that the second H-like 
iron component is created in a warm (T ~ 10 6 K) gas within a few 
100s R s of the central source. 

Finally to check our results we fitted the the ASCA spec- 
trum above 1 keV (SIS0 and GIS2, sequence number 71028030, 
from the Tartarus database) with an ionized disc reflection model. 
The fit was good (\ 2 /do f =399. 3/378) and the ionization parame- 
ter was consistent with that observed in the BeppoSAX spectrum, 
log (,asca = 1.08±o;oi- A narrow line component also is really 
marginally required (Pf=92 per cent) to fit the Fe line. The flux 
of this component was the same in the two observations, while the 
flux in 2-10 keV was higher by 30 per cent during the SAX long 
look. Nevertheless the uncertainty on the line intensity rends any 
result about variability inconclusive. 
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NGC 7469. Ionized disc reflection model 




2.01 
Photon Index 



Figure 6. 68, 95 and 99 per cent confidence contours Reflection fraction vs 
photon index for the ionized disc plus gaussian line model. The "+" denotes 
the best fit values of R and T. 
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